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PHASE MEASURING METHOD AND APPARATUS FOR MULTI- FREQUENCY INTERFEROMETRY 

The present invention relates generally to phase ' 
measuring sensors for measuring the phase of 
5 electromagnetic waves and, in particular, though not 
exclusively, to methods for determining absolute fringe 
order in fringe patterns obtained in such sensors. 

The phase of electromagnetic waves projected onto an 
object can be affected by various measurands of the object 

10 such as deformation [1], vibration [2], and refractive 
index variation due to density changes. Phase measuring 
sensors are often used for singe point or whole field 
prof ilometry, for example profiling of three-dimensional 
objects. There exist a host of interf erometric based 

15 techniques for single point and whole field prof ilometry, 
from synthetic aperture radar (SAR) [3] to fringe 
projection [4]. In fringe projection techniques, a 
projected wavelength is generated by projection of a known 
fringe pattern (sometimes referred to as a ^^fringe map'') 

20 onto the object of investigation, under a certain angle of 
incidence. In interf erometric based sensors the third 
dimension of an object is coded in an interf erogram or (in 
fringe projection techniques) in a deformed fringe pattern. 
In the following it will be understood that a deformed 

25 fringe pattern can be analysed, and phase information 

extracted therefrom, in the same way as an interf erogram. 

In an interference pattern there are two elements that 
contribute to the phase measurement dynamic range: the sub- 
fringe phase resolution and the number of fringes (fringe 

30 orders) spanned by the measurement. There is typically a 
simple function relating the measured optical fringe phase 
to the desired measurand. The process of phase measuring 
produces sub-fringe resolution, typically 1/100^*" to 1/1000^^ 
of a fringe. The sub-fringe resolution is calculated via 



wo 2004/020937 




PCT/GB2003/003744 



2 



either phase stepping [5] or Fourier transform techniques 
[6] . However, the interf erometric phase is calculated using 
an inverse trigonometry function with principal values over 
the range to +7C at best. Therefore the required phase 
5 information is ^wrapped' into an interval with sharp 
discontinuities in the data at the edges of that interval 
which must be spatially or temporally unwrapped to obtain 
the fringe order information. This problem is typically 
resolved by phase "'unwrapping'' (alternatively known as 

10 fringe counting) . Conventional interf erometric analysis for 
single point data relies on a temporal scan to measure the 
relative phase change in going from one state to another. 
In the case of whole field data, a spatial unwrapping is 
achieved using an appropriate algorithm giving relative 

15 information on the state (or change of state) across the 
image field [7] . 

A big disadvantage of conventional phase measuring 
systems is that only '^relative'' information can be 
obtained. Many industrial measurements require absolute 

20 information, for example in the measurement of range. In 
other cases absolute data is needed to overcome 
implementation issues which mean it is not possible to 
fringe count sufficiently rapidly to track from one state 
to another and so an absolute measurement of each state is 

25 necessary. In whole field sensors, a generic problem common 
to most interf erometric techniques is the determination of 
absolute fringe order in interf erograms/deformed fringe 
patterns containing phase discontinuities or spatially or 
temporally discrete samples. In such cases it is difficult 

30 to unambiguously determine fringe orders - giving absolute 
data - in the spatial unwrapping process. The spatial 
unwrapping process relies on the field being spatially 
contiguous i.e. it cannot be applied to determine relative 
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phases for spatially separate objects or objects containing 
a discontinuity. 

One approach to the unambiguous determination of 
fringe orders has been through the use of two wavelength 
5 interf erometry techniques [see H Zhao, W Chen, Y Tan, Phase 
unwrapping algorithm for the measurement of 3D object 
shapes. Applied Optics, 33, 4497-4500 (1994)]. Here two 
wavelengths are heterodyned to generate a beat wavelength. 
In two wavelength interf erometry using wavelengths and 

10 ^1, with Xo < Xx, the unambiguous measurement range is 
increased (compared with single wavelength interf rometry) 
to the synthetic wavelength Aoi at the beat frequency, where 
Aoi = XqXi/ {Xi-Xo) . From this, absolute data can be extracted. 
The ability to calculate a fringe order from measurements 

15 at two wavelengths removes the need to spatially unwrap, 
assuming that the heterodyne process is robust. However, 
there is a finite, non-zero, phase measurement noise in two 
wavelength heterodyne systems and as this phase noise 
increases the number of fringes which may be successfully 

20 heterodyned decreases. This means that two wavelength 
heterodyne sensors tend not to be very robust, in terms of 
giving unambiguous absolute data. 

Another well established research area for absolute 
phase measurement is fringe projection for profilometry 

25 where the projected wavelengths produced can be varied 
relatively easily by varying the projected fringe pattern 
or the angle of incidence. Therefore techniques based on 
temporal phase unwrapping may be applied. Huntley et al 
introduced a temporal phase unwrapping approach using a 

30 reverse exponential series of projected fringes given by: 
5, 5 -1,5 -2, 5- -4,..., 5/2 [8] . This is the subject of International 
Patent Application Publication No. WO 97/36144. This 
technique allows 2^"''" fringes to be counted absolutely where 
N is the number of fringe frequencies used. Here, the 
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unwrapping is performed between consecutive phase 
measurements to scale the fringe order calculated at the 
subsequent wavelength. Absolute fringe order is obtained as 
the unwrapping is performed at each pixel independently 
5 along the time axis. The number of projected fringe 
frequencies required is (logj(j)+l). For each projected fringe 
frequency four phase stepped images are obtained to 
determine the wrapped phase values in the interval -n to h-tc. 
For example, to measure over a range of 32 fringes 6 sets 

10 of 4 frames are needed, and for 128 fringes 8 sets of 4 
frames. Therefore, considerable time is required to obtain 
the image frames, and as each image may contain >1MB of 
information a significant data processing problem is 
generated. Furthermore, a mechanism must exist to generate 

15 the wide range of wavelengths required- This would be a 
particular problem if a technique other than fringe 
projection were being used as then it would be the optical 
wavelength which must be varied (not simply a synthetic 
wavelength) . 

20 It is an object of the present invention to avoid or 

minimise one or more of the foregoing disadvantages. 

According to a first aspect of the invention there is 
provided a multi-frequency interf erometry method for 
measuring the absolute phase of an electromagnetic wave, 

25 wherein the method comprises: 

selecting a target measurement range, L, within which 
absolute phase measurements are desired to be made, said 
phase measurements being related to a desired measurand of 
an object; 

30 determining a level of phase noise which will be 

present in wrapped phase measurements which will be made; 

for the selected target measurement range, and the 
determined phase noise level, calculating an optimum 
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number^ of wavelengths of electromagnetic radiation to 

be used in the multi-frequency interf eroihetry, where N>3; 

selecting an optimum series of values of said N 
wavelengths to achieve optimum noise immunity in wrapped 
5 phase measurements to be made; and 

carrying out multi-frequency interf erometry at the 
selected values of the N wavelengths so as to make at least 
one wrapped phase measurement at each of the N wavelengths, 
which wrapped phase measurements are processed to obtain an 
10 absolute phase measurement related to the desired measurand 
of the object. 

The processing of the wrapped phase measurements may 
include one or more of heterodyne processing, Fourier 
series processing, recursive unwrapping and iterative 
15 unwrapping. 

The optimum number, N, of wavelengths is preferably 
the minimum number of wavelengths required to obtain 
unambiguous phase measurements in the selected target 
measurement range, for the determined phase noise level. 

20 The method may include proposing a measure of process 
reliability associated with the determined phase noise 
level, corresponding to a known probability of success in 
fringe order identification. For example, a process 
reliability of 6a may be proposed, where a is the standard 

25 deviation noise in a function (for example in the form of a 
fringe pattern) to be analysed: this function could, for 
example, be a discrete level heterodyne function generated 
from the wrapped phase measurements made at two different 
wavelengths (this discrete level heterodyne function will 

30 be used to obtain absolute phase data) . The proposed 
process reliability is preferably used to calculate the 
minimum number, N, of wavelengths required to obtain 
unambiguous phase measurements in the selected target 
measurement range • 
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The selected values of said N wavelengths may define a 
geometric series or, alternatively, may be combined in 
post-processing techniques, for example heterodyne 
processing, to generate a geometric series of synthetic 
5 wavelengths. Where heterodyne processing is used, 
preferably the selected values of the N wavelengths are 
chosen such that only one heterodyne operation is required 
to generate each synthetic wavelength in the geometric 
series. If more than one heterodyne operation is required 
10 to generate a desired synthetic wavelength it will be 
appreciated that measurement noise levels will increase. In 
a preferred embodiment, the values of said N wavelengths 
may be selected in accordance with the following equation: 

15 lAx = lAo - (l/>uo)^"^^'''^(l/L)^-^^''-^ (Equation A) 

where i=l,...,N, where N is the number of wavelengths, A.o is 
the wavelength associated with the largest frequency, and Xi 
is the wavelength of the ith frequency^ and L is the target 

20 measurement range. However, there are also other possible 
series of projected wavelengths which will generate the 
desired geometric series of synthetic wavelengths. 

The inventive method may be applied in many different 
types of interf erometry, for example, in single point and 

25 full field prof ilometry, synthetic-aperture radar (SAR) 
interferometry, multi-aperture synthesis techniques, and 
fringe projection techniques. For example, the method may 
be used for measuring the three-dimensional shape of an 
object which may have surface discontinuities. In this 

30 case, the target measurement range is preferably the range 
over which (absolute) depth measurements relating to the 
profile of the object are desired. (Alternatively, the 
target measurement range may be the desired field of view 
of the imaging system in which the method is applied.) The 
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Step of carrying out the multi-frequency interf erometry may 
conveniently comprise: 

recording a series of fringe patterns obtained when 
the object is illuminated with electromagnetic radiation at 
5 each of the selected N optimum wavelengths; 

processing the recorded fringe patterns so as to 
obtain wrapped phase measurements in the form of a wrapped 
phase map for each of the N optimum wavelengths, and 
processing the wrapped phase maps to determine absolute 

10 fringe orders in the fringe patterns. The absolute fringe 
orders may be calculated using heterodyne processing and/or 
iterative unwrapping of the wrapped phase maps. The 
determined absolute fringe orders can then be used to 
unambiguously unwrap the wrapped phase maps, from which 

15 unwrapped phase maps a three-dimensional profile of the 
object can be compiled. 

It will be appreciated that, depending on the 
particular interf erometry technique being used, there will 
be various ways of illuminating the object with radiation 

20 at the desired optimum wavelengths. For example, in fringe 
projection based interf erometry, white light illumination 
is used to project a predetermined fringe pattern on to the 
object to be measured. In this case, it will be appreciated 
that the object is in effect illuminated at a synthetic 

25 wavelength (hereinafter referred to as the ^''projected 

wavelength"), this projected wavelength being determined by 
the chosen fringe pattern projected onto the object. In the 
present invention, if the multi-frequency interf erometry is 
carried out using fringe projection, the selected series of 

30 N optimum wavelengths are therefore different projected 
wavelengths at which the object must be illuminated. In 
other forms of interf erometry the selected series of N 
optimum wavelengths may be actual different wavelengths of 
light, normally in the form of laser beams, with which the 
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object must be illuminated (usually sequentially) . Thus, 
the recorded fringe patterns may be interf erograms 
(produced in phase measuring sensors, including 
interferometers and single point and full-field systems) or 
5 may be deformed fringe patterns produced in fringe 
pro j ection interf erometry . 

If the multi-frequency interf erometry is carried out 
using fringe projection, the optimum values of the 
projected wavelengths may be selected in accordance with 
10 the following re-written version of Equation A: 

N^^ = N^^ - (a^^^ , for / = 1,2,.... A ( Equation B ) 



where X is the number of projected wavelengths, N^^ is the 
number of projected fringes in the largest fringe set, and 
15 is the number of projected fringes in the fringe set. 

Alternatively, where the number, N, of wavelengths is 
three, the selected values of the three wavelengths may be 

Preferably, at each said optimum wavelength a series 
20 of fringe patterns is recorded. Most preferably, a series 
of phase-stepped fringe patterns is recorded at each of the 
N wavelengths with which the object is illuminated. The 
phase-stepped fringe patterns can be used to create a 
wrapped phase map for each said optimum wavelength. 
25 Alternatively, the wrapped phase maps may be obtained using 
Fourier transform processing of a single fringe pattern 
recorded at each said optimum wavelength. 

The object may be illuminated separately and 
sequentially with each of the selected optimum (real or 
30 projected) wavelengths of electromagnetic radiation. 
Alternatively, the object may be illuminated with white 
light and a plurality of the desired fringe patterns may be 
captured simultaneously by recording them with an image 
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detector which can simultaneously capture image data at the 
selected optimum wavelengths. For example, if the selected 
optimum wavelengths lie in discrete blue, green and red 
colour bands, a colour camera can be used to separately but 
5 simultaneously capture image data at red, green and blue 
frequencies* In another possible embodiment, the inventive 
method may be used in high speed single point ranging 
applications- In this case the desired measurand is of 
course absolute range. In one possibility, the multi- 

10 frequency interf erometry may be carried out by illuminating 
the object using a broadband femtosecond laser. Preferably, 
the output of the laser is split into two parts: one part 
illuminates the object and the other part is used as a 
reference arm. The interference of the object and reference 

15 light produces an interf erogram in the from of a time- 
varying intensity signal which may conveniently be recorded 
using a single point detector, for example a photodiode. 
The spectral range of the broadband pulses emitted by the 
laser is chosen so as to include the desired optimum series 

20 of values of said N wavelengths, for the target measurement 
range, L, selected in accordance with the inventive method. 
Other unwanted wavelengths can be filtered out, so that 
only interferograms at each desired optimum wavelength are 
recorded. Phase measuring means may be provided to enable a 

25 wrapped phase measurement to be obtained for every, or 
nearly every, femtosecond pulse emitted by the laser. The 
absolute fringe orders in the interferograms may be 
determined by, for example, heterodyne processing or 
iterative unwrapping of the wrapped phase measurements. 

30 From the above it will be generally appreciated that 

the invention provides a novel strategy for absolute fringe 
order identification in multi-wavelength interf erometry 
based on optimum selection of the wavelengths to be used. 
As will be described in detail below with reference to the 
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preferred embodiment a theoretical model of the process 
has been developed which allows the process reliability to 
be quantified. The methodology produces a wavelength 
selection which is optimum with respect to the minimum 
5 number of wavelengths required to achieve a target dynamic 
range. Conversely, the maximum dynamic range is produced 
from a given number of optimally selected wavelengths 
utilised in a sensor. Thus, according to a further aspect 
of the invention there is provided a method of measuring 

10 absolute fringe order in a phase measuring sensor, wherein 
the method comprises: 

selecting a number, N, of wavelengths of 
electromagnetic radiation to be used to illuminate an 
object, where N>3; 

15 determining a level of phase noise which will be 

present in wrapped phase measurements which will be made; 

for the selected number, N, of wavelengths, and the 
determined phase noise level, selecting optimum values of 
the N wavelengths so as to achieve a maximum measurement 

20 range in which wrapped phase measurements can be made 
relating to a desired measurand of the object to be 
illuminated; 

recording a series of fringe patterns obtained when 
the object is illuminated with electromagnetic radiation at 
25 each of the selected N optimum wavelengths; and 

processing the recorded fringe patterns to obtain at 
least one wrapped phase measurement at each of the selected 
N optimum wavelengths, and processing the wrapped phase 
measurements to determine absolute fringe orders in the 
30 fringe patterns. 

According to another aspect of the invention there is 
provided a multi-frequency interferometer apparatus for 
shape measurement, comprising: 
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fringe projection means for generating a projected 
wavelength of illumination by projecting a predetermined 
fringe pattern on to an object; and 

image capture and recording means for capturing and 
5 recording a deformed fringe pattern obtained when said 
predetermined fringe pattern is projected on to the object, 
the image capture means being disposed at an angle to the 
direction of illumination of the object with the projected 
fringe pattern; and 

10 data processing means for processing recorded deformed 

fringe patterns so as to obtain phase measurements 
therefrom; wherein 

the fringe projection apparatus is variable such that 
an optimum series of projected wavelengths may be 

15 generated, the values of the wavelengths being such that, 
for a known level of phase noise in the apparatus, and a 
chosen number of projected wavelengths in said series, 
absolute fringe orders in the deformed fringe patterns are 
measurable over a maximum measurement range. 

20 The fringe projection means may be any known means for 

producing projection fringes. For example, the fringe 
projection means may comprise a coherent fibre fringe 
projector for producing a pattern of Young's fringes across 
the object, or may comprise a spatial light modulator based 

25 fringe projector, or imaging means for imaging grating 
patterns on to the object. 

It will be appreciated that the new concept introduced 
for optimum wavelength selection is readily scalable, i.e. 
from a three wavelength system to a four wavelength system, 

30 from four wavelengths to five etc. Therefore, whilst the 
previous technology described in W097/3614 4 and reference 
[8] allowed 2^"^ fringes to be counted absolutely where N is 
the number of fringe frequencies used, the new technique 
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scales as n**"-^, where n is an arbitrary real number limited 
by phase noise of a practical interferometer. 

Preferred embodiments of the invention will now be 
described by way of example only and with reference to the 
5 accompanying drawings in which: 

Fig.l illustrates graphically the three discrete level 
heterodyne functions for an optimum A-X process; 

Fig. 2 illustrates graphically the three discrete level 
heterodyne functions for an optimum process with a 

10 total of 80.765 fringes; 

Fig. 3 is a graph of Fringe Order against Sample Number 
across an area detector^ where the Sample Number is a 
function of position across an object, calculated for the 
same optimum process as Fig. 2, having a total of 80.7 65 

15 fringes; 

Fig. 4 (a) is a schematic diagram of an experimental 
set-up for fringe projection shape measurement; 

Fig. 4(b) is a schematic diagram illustrating a typical 
fibre fringe projector for use in the experimental set-up 
20 of Fig. 4 (a) ; 

Fig. 5 is an experimentally obtained fringe order map 
for optimum heterodyne processing; and 

Fig. 6 is a schematic diagram of a high speed single 
point ranging sensor. 

25 

A detailed analysis of the theoretical basis of the 
invention will now be described, followed by a description 
of various practical embodiments. This analysis is 
presented with reference to fringe projection 
30 interf erometry, but it will be appreciated that the 
analysis is equally applicable to all forms of 
interf erometry in which it is desired to determine absolute 
fringe order in an interf erogram. 
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Theoret:ical Developmezi't 

Two Wavelength Interferometry 

In two wavelength heterodyne fringe projection 
interferometry to eliminate step height ambiguity the 
5 difference in the number of fringes projected across the 
field of view must be <1. Let the number of projected 
fringes across the field be A^^ for the wavelength; then 
Nj^^-N^^<\ with Nj^,>Nj^^. The difference in wrapped phases 
calculated at the two wavelengths, a heterodyne function, 

10 can be expressed as a phase within the interval -n to n and 
will consist of a monotonic ramp across the image. A 
convenient representation is to calculate a discrete phase 
level for each of the fringe orders corresponding to NJ^^ . 
This is obtained by subtracting a scaled version of the 

15 wrapped phase at iV^, from the heterodyne function - the 
scaling factor being given by (a^^, - A^«)/^ii • practice, such 

a discrete level heterodyne function can only identify a 
limited number of fringe orders owing to the presence of 
phase noise* Each //^ is known, as it is set in a white 

20 light system or can be measured in coherent systems [9] . 

Each wrapped phase measurement contains phase noise, which 
is modelled as a Gaussian distribution with zero mean and a 
standard deviation of cr^ [10] . The heterodyne function then 

contains noise with a standard deviation given by V2cr^ . We 
25 define a process robustness of 6a , corresponding to a 
probability of 99.73% success in fringe order 
identification within the measurement system. Therefore, 
the discrete phase levels must be separated in phase by at 
least 6yf2a^ . Hence, for 6a reliability the number of fringes 
30 which can be correctly identified is limited by: 



wo 2004/020937 PCT/GB2003/003744 

14 

Defining the measurement dynamic range as the product of 
the phase resolution and the number of fringes successfully 
numbered. Therefore, for a phase resolution of l/IOO//rof a 
fringe, 12 fringes can be numbered to 6a reliability, 
5 giving a dynamic range of 1200, which is insufficient for 
most engineering applications. 

Optimally Selected Wavelengths in N Wavelength Interferometry 

The introduction of a third projected wavelength allows the 
10 generation of two independent heterodyne functions 

containing N^^, and N^^^ discrete levels, where N^^xN^^=Nj^^ . 

The discrete phase levels in each modified heterodyne 

function will be separated by: 

3^,=27r/N^, , and ^^^=2;r/N^, . (l) 
15 In this arrangement it is found that as increases 

decreases. Therefore, for maximum overall reliability in 

fringe order identification the optimum set of projected 

fringe wavelengths is given by the symmetrical arrangement 

where ^ou=^di2 ^rid ^du^^du^^^x * the number of discrete 

20 levels is not balanced the effect is to increase one of iV^, 
and Ar^2/ thereby bringing ,9^^, or ,9^„ closer to the noise 
limit. Hence, from equation (l), with an optimised three 
wavelength approach the condition for 6cr reliability fringe 
order identification is given by. 

The number of fringes which may be heterodyned reliably in 
an optimum three wavelength heterodyne set-up is therefore 
the square of that for the two wavelength case for the same 
phase measurement noise. 
30 Equations 1,2 are used to define general expressions 
relating the numbers of projected fringes at the N 
wavelengths. By assuming N^^>Nj. and 
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with ^oii =-^ii,.2 =- = A^£«/ / ^or optimum noise immunity, it is here 
proposed that the numbers of projected fringes should be 
selected in accordance with the following expression: 

5 for I = 1,2,..., A (3) 

where: A number of wavelengths, Ny,^ number of fringes in the 
largest fringe set, N^^ number of fringes in the /'* fringe 
set. The term ^,,=0 is included in order to generalize the 
expressions that follow. The general formulas to calculate 
10 the fringe order are given by: 

^^'^^'^^^ N -N ^-^0... (4) 



2^ 



(5) 



for 1 = 1 ,>l-l , where DL^ is the /'* discrete level function, 

IDL. is the i'* integer discrete level, H^^ is the heterodyne 

15 between the 0 and i"* wrapped phase map and NINT denotes 
taking the nearest integer. The recursive relationship for 
the integer discrete level functions is initialised by 
setting JDLq=0. The fringe order for the wrapped phase map 
with N^^ fringes is given by/DZ^ from equation (s) . 

20 For example using four synthetic (i.e. projected) 
wavelengths with N^^=64, then iV,, =63, n^^=60 , AT^, =48, Fig. 1 

shows the discrete level functions resulting from a 
simulation of the process. The phase difference between 
discrete levels (y-axis) is equal in each plot as expected 

25 for an optimum wavelength configuration. In this case Nxi = 
64 and this means the discrete level functions contain 1 
fringe with 4 discrete levels on it (1*4=4 fringe orders), 
4 fringes each with 4 discrete levels on them (4*4=16 
fringe orders), and 16 fringes each with 4 discrete levels 

30 on them (4*16=64 fringe orders) . So we have the pattern 1, 
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4, 16, 64, or to be generic: (Nu )^^^, (Nxi )^^^, Nxi , where 
we note that we are dealing with cube roots which is the 
number of wavelengths-1 . The number sequence here defines a 
geometric series, where the factor relating neighbouring 
5 terms is Ndl- The geometric series can be expressed as 

(^/o)^' for / = l,2,...,A, where A is number of wavelengths. 

A multiplicative intensity noise of 2.5% has been 
applied to the cosine intensity fringes which is 
representative of scientific CCD cameras and is a primary 
10 source of phase noise. 

It is here proposed that the number of fringes that may be 
numbered to 6cr reliability for a system with A wavelengths 
is given by (from equation (2)above ) : 



15 The reliability of the process only depends on the first 
and last wrapped phase maps with all intermediate values 
for other wavelengths cancelling out (see equations 4 and 
5) . As the last wavelength corresponds to zero fringes, 
rather than performing a phase measurement here it is more 

20 accurate to set the wrapped phases for this wavelength to 
zero. Hence the process reliability depends only on the 
error in the wrapped phase map at Nj^^ . Therefore, as the 

number of wavelengths used is increased the process 
reliability remains unchanged and the number of fringes 
25 that may be labelled increases as given by equation 6. 

An alternative representation of equation 6 is useful 
for some applications and obtained by noting that i//L^,. =iVyy 

where: is the wavelength associated with the i^^ 

frequency and L is the desired unambiguous measurement 
30 range. Equation 6 would then be re-expressed as: 




(6) 



wo 2004/020937 




PCT/GB2003/003744 



17 



Fig.l shows simulated (i.e. computer generated) phases 
calculated at the synthetic (i.e. projected) wavelengths 
for optimum four-frequency interf erometry with L/X^-M . 

In any measurement situation the phase noise will be 
5 measured and hence known, and the target measurement range, 
L, will be known. Then equation 6 can be used to determine 
the number of wavelengths needed to make ^xx large enough to 
span the desired measurement range, in that: 

Nxi = L/A,o / where Xq is the wavelength associated with the 
10 largest frequency. The result of this will be the minimum 
number of wavelengths, X, needed to achieve the target 
measurement range for the particular sensor system. 

It should be noted that in the example of Fig.l, the 
value (A^yj,)*"' is an integer and therefore the discrete levels 
15 are equal in phase values for each discrete level function. 
For any real positive number, {j^fy~\ the algorithms function 
correctly. For example, AT^^ =80.765 and using four projected 
fringe wavelengths the discrete level functions are given 
in Fig. 2 (again showing simulated phases). The recursive 
20 unwrapping process correctly unwraps successive beat 
fringes to obtain the fringe order, as seen in Fig. 3. 

For a typical four wavelength system cr^ =2;r/85 and the 

number of fringes that can be correctly identified with 60- 
reliability is >1000, giving a dynamic range of >85000. The 
25 following table (Table 1) shows values for N^^ and dynamic 
range for a range of X and . 



Phase 


Optimum 


Optimum 


Noise 


3-Z Heterodyne 


4'A Heterodyne 


(rms. 


Nj.^ Dynamic 


N^^ Dynamic Range 


radians) 


Range 




InjSO 


35 1,736 


205 10,230 
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2ff/100 


139 


13,889 


1,637 


163,700 




556 


111,111 


13,095 


2,619,000 




2,222 


888,889 


104,707 


41,902,800 



Table 1: Numbers of fringes and dynamic range obtained with 
99.73% reliability for optimum 3-X and heterodyne 
processing. 



5 It will be understood that although in the above 
description we have proposed 6a process reliability, the 
process reliability may be defined by other values as 
desired, for example 8a reliability, in which case equation 
(6) becomes: 

Alternative Strategies for Generating Optimum Wavelengths 

There are alternative ways to calculate the numbers of 
fringes used at each wavelength that result in optimum 
15 configurations in the discrete level functions. For example 
using the three projected wavelengths given by Nj^, N^^ - ^N^^ , 

yjNj.^ -1 , these can be heterodyned together as follows to 

obtain the desired geometric series 1, ^Z^/, / : heterodyne 

N^^ with N^^-^jW^ to get '^N^^, heterodyne this with ^^-1 to 

20 get 1. Equations 4 and 5 may than be applied as before. It 
is expected that the noise floor is worse in this case 
(than for wavelengths selected in accordance with Equation 
3) as it relies on two heterodyne operations to get one of 
the synthetic wavelengths (the 1) in the geometric series . 

25 There are also many other possibilities to arrive at the 
desired geometric series 1, ^^/^f the only difference 

in the various ways is the number of heterodyne operations 
required in each case. However, the noise level will 
increase the greater the number of heterodyne operations 
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that are required to obtain any synthetic wavelength in the 
desired geometric series. 

In fact, if the maximum number of heterodyne 
operations to get to a member of the geometric series is 
5 denoted by r, the measurement range equation 6 can be re- 
expressed as: 

For example, in the three wavelength case, if the original 
projected wavelength selection is 100, 99, 90 a single 
10 heterodyne is needed to get to each term in the desired 
geometric series 1,10,100 and hence r=l, leading to 
equation 6. 

In fact, in the case of fringe projection using three 
wavelengths, one could capture directly a fringe pattern 

15 with the projected wavelengths corresponding to 1, Va^^^, and 
Nj.^ fringes and apply iterative unwrapping to these 
directly. This would be instead of the more usual capture 
of N^^ , N^^-^N^^ (as proposed by equation 3), then 

heterodyne (i . e . ''^beat") A^,^ with to get 1, and beat 

20 Nj.^ with Nj.^-^Nj.^ to get the ^N^^ and then apply iterative 
unwrapping to the derived data at 1, ViV^^ , and N^^ fringes. 

Capturing directly a fringe pattern at the series 1, 
^Nj,^, and N^^ fringes thus corresponds to the case r=0 in 

equation 7 above. In this case it will be appreciated that 
25 discrete level functions can be formed directly in the 
processing stage without heterodyning first. The discrete 
level function is obtained by subtracting a scaled version 
of the wrapped phase at the next higher frequency in the 
series, e.g. for 3 wavelengths to get ^N^^ levels on the 

30 single unambiguous fringe the phase measured with '^Nj.^ 

fringes is multiplied by and subtracted from the 
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phase obtained with a single fringe. The discrete level 
functions may then be iteratively unwrapped. The geometric 
series 1, ViV^^ , N^^ still provides an equal factor of ^N^^ in 

going from 1 to ^N^^ to N^^ . As these factors are equal, 

5 when the discrete level functions are calculated they will 
contain equal numbers of discrete levels per fringe and 
hence the process reliability is equal in the discrete 
level functions and the reliability of the process overall 
is optimised. From this it will be understood that in some 
10 possible embodiments of the invention the projected 
wavelengths may be chosen so as to match exactly the 
desired geometric series i.e. the series defined by 

O^/o)^' for f = l,2,...,A, where A is number of wavelengths. 

15 Modifications to the Theoretical Description 

To allow for experimental error in obtaining exactly the 
number of fringes required across the measurement range, 
the wavelengths would be selected close to their 
theoretical values but in such a way as to define an 

20 unambiguous measurement range slightly larger than that 
required. For example, the second wavelength (1=2, in 
equation 3) may be modified by up to 30% in order to 
guarantee forming less than one beat fringe across the 
desired measurement range. Therefore, in practical systems 

25 it would be common to use Nj,^ , A^,, -0.7, Nj.^ - ^N^^ for a three 
wavelength system. 

Experimental Demons 'bratlon 

A whole field, triangulation based shape measurement 
30 system is given as the exemplar for the above theoretical 
analysis. In this system a coherent fibre fringe projector 
10 is used to produce a pattern of Young's fringes 8 across 
a test object 12. A CCD camera 14 is then placed at an 
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angle 9 to the illumination direction to view the fringes^ 
as shown in Fig. 4(a). The CCD camera 14 is connected to a 
computer system 16 for implementing the data processing 
scheme (in particular, the heterodyne operations and the 
5 recursive unwrapping algorithms specified by equations (4) 
and (5) above) and control of the data acquisition process. 
The computer system includes a memory for recording fringe 
patterns captured by the CCD camera 14. The computer system 
is programmed to carry out the phase unwrapping and fringe 

10 order identification processes as herein described in 
accordance with the invention. Increasing the number of 
projected fringes or the angular separation 9 between the 
CCD camera 14 and the fringe projector 10 increases the 
sensitivity to the depth (z) of the test object 12. 

15 Fig. 4(b) illustrates in more detail the features of the 
fibre fringe projector 10. Like parts in Figs. 4 (a) and(b) 
are referenced by like numerals- In known manner, this 
comprises two output fibres 20,21 having their ends placed 
side-by-side with a spacing therebetween. One of the output 

20 fibres 20 was wrapped around a cylindrical PZT 24 to allow 
the optical phase to be modulated. A servo control system 
26 allows the phase of the projected fringes to be 
stabilised by monitoring the reflections from the distal 
ends of the output fibres at the 4*^^ arm of a directional 

25 coupler 30. The servo also enables accurate 90° phase steps 
to be obtained [11] - A second CCD camera 32 was 

incorporated into the fringe projector to measure the fibre 
separation, and hence the number of projected fringes. The 
fringes were sampled directly onto the second CCD camera 32 

30 via a polarising beamsplitter 34. A measurement of the 
phase distribution across this CCD camera 32 gives the 
fibre separation and hence the number of projected fringes 
[9] . The fibre separation can be varied by means of a 
linear traverse (not shown) on which one of the fibre ends 
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is mounted. We have demonstrated a resolution of <10 nm in 
the fibre separation measurement compared to a 50 nm 
requirement predicted by a simulation for three projected 
wavelengths. To demonstrate the process, the flat side of an 
5 object was assessed initially as this gave easy 
verification of the fringe order calculation. Three 
wavelengths were used (from equation 3 above) with 100/ 99 
and 90 projected fringes. Fig. 5 shows the result of 
applying equations 4 and 5 to the three wrapped phase maps 

10 obtained. The central fringe is automatically identified by 
the heterodyne process and is coloured black. The fringe 
orders are clearly identified as a repeating scheme of 6 
colours either side of the central fringe. 

The phase resolution obtained was estimated to be 

15 1/80^^ fringe principally limited by the random speckle 
phase which contributes to the interference phase. With 
o-^=2;r/80 and ;\^,„=100 equation 2 can be manipulated to 

determine the fringe order numbering reliability expected 
from theory which evaluates to 99.53%. A local 

20 neighbourhood check applied to the central flat area of the 
experimental data in Fig. 5 shows the number of pixels 
giving the correct fringe order to be 99.52%. Therefore, a 
validation of the process and the underlying theory has 
been demonstrated . 

25 The new process is equally applicable to white light 

fringe projectors based on spatial light modulators or on 
custom gratings produced on glass substrates. Compared to 
the reverse exponential temporal unwrapping scheme [8], to 
measure over --100 fringes requires 8 wrapped phase maps to 

30 be obtained. Therefore, the use of the new algorithm 
enables a reduction in data acquisition time and data space 
of >60%. The dynamic range obtained in this measurement was 
1:8000. 
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The above described embodiment provides a robust 
measurement process based on multiple wavelengths, phase 
measurement and heterodyne processing. It is shown that 
given a certain phase measurement noise it is possible to 
5 maximise the number of fringe orders unambiguously 
identifiable and therefore the measurement dynamic range. 
Each fringe order can be identified as an absolute number. 
A robustness measure based on Gaussian statistics is 
derived for the process. 

10 It will be appreciated that the inventive process is 

applicable to any form of interf erometric sensor where 
phase measurements are available at defined sensitivities 
(wavelengths), not simply only to fringe projection based 
phase measuring sensors. For interf erometric sensors not 

15 based on fringe projection, the equivalent of the optimum 
noise immunity condition = AT^^ =- = -^/>/^ proposed above for 

fringe projection sensors is: 

A01/A02 = A02/A03 = = Ao, (8) 

where Aoi, A02, A03,..., Ao, n-i are the synthetic (beat) 
20 wavelengths formed when Xq is heterodyned with A,i, X2, Xn-i 
respectively, e.g. Aoi is the phase difference between phase 
measurements at A.o and X-i and is defined by: Aoi = X.oXi/(A-i- 
X.o) • In this case the equation (3) is re-expressed as: 

25 l/>.i = lAo - (lAo)^"^^''"Ml/L)^"^/'^-^ (Equation A) 

where i=l,...,N, where N is the number of wavelengths, Xq is 
the wavelength associated with the largest frequency, and Xi 
is the wavelength of the ith frequency, and L is the target 
30 measurement range. However, again it will be appreciated 
that alternative series of wavelength values are possible 
which can be combined in post-processing to produce the 
desired geometric series of synthetic wavelengths 
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satisfying equation 8 above. With reference to equation (7) 
above, the r=0 case is also possible (i.e. the wavelengths 
are chosen to match the desired geometric series, without 
processing heterodyne operations being required ) as long 
5 as the required range of wavelengths is available from the 
wavelength sources being used and are detectable using 
sensors. In practice, the visible wavelength range will not 
be broad enough to offer an r=0 solution, but an example 
where it would be possible is in radar systems where the 
10 wavelength range available is orders of magnitude broader. 

Examples of alternative inventive embodiments (not 
based on fringe projection) are: 

a) High speed, single point ranging. In this application 

15 the aim is to monitor one point on a target object at 

high speed . Short pulse, pico and femtosecond, lasers 
produce a pulse of light over a broad range of 
wavelengths (as the pulse length shortens the 
wavelength range increases) . Fig. 6 illustrates an 

20 example embodiment. The light from a pico or 

femtosecond laser 40 is split into two parts (using a 
first directional coupler 41) , one part goes to the 
object 50 and the other acts as a reference beam 42 in 
a variable delay line 44 (for path length matching) . 

25 Light scattered by the object 50 is then combined with 

the reference light (using a second directional 
coupler 43) • For each wavelength required, a filter 46 
is used to isolate that wavelength from the broad 
bandwidth pulse, with the filtered light detected on a 

30 single point detector 48 (a suitable example would be a 

photodiode) . The output from the detector 48 is passed 
to a high speed data logger (not shown) . In practise, 
a pair of detectors would be used in quadrature for 
each wavelength channel in order to calculate the 
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optical phase. Thereby a phase measurement is obtained 
on each pulse from the laser. So the ^interf erograms' 
(equivalent in purpose to the deformed fringe patterns 
used in the above-described fringe projection 
5 profilometry embodiment of Fig. 4) are mapped out in 

time from each detector as an intensity signal. The 
absolute phase information is obtained by post- 
processing the interf erograms obtained at each 
wavelength in a similar manner to the post-processing 
10 of the wrapped phase maps obtained in the Fig. 4 

embodiment, for example using heterodyne processing to 
generate data at the desired optimum geometric series 
of synthetic wavelengths, to extract the absolute 
phase data (from which absolute range measurements are 
15 calculable) . 

b) Colour camera application. In this embodiment 
simultaneous illumination of the object occurs at three 
wavelengths. The three wavelengths are chosen to reside in 
discrete colour bands, e.g. the red, green and blue 
20 components of a colour camera (having a known arrangement 
of red, green and blue pixels) with zero or minimal leakage 
of one component into the neighbouring colour band A colour 
camera provides simultaneous acquisition of the images in 
each of the colour bands. The colours may simply be used as 
25 a means of separating out the multi-wavelength image 

obtained at the colour camera in which case it would be 
preferable to obtain directly the data at the 1, Vn, and N 
fringes, where N is the number of fringes in the highest 
resolution fringe set. Alternatively, the wavelengths would 
30 be selected to give (with some processing e.g. heterodyne 
operations) the geometric series defined by 1, Vn, and N 
fringes. Phase data captured at these wavelengths could be 
heterodyned to get to the 1, Vn, and N fringes and hence 
iteratively unwrapped. This latter approach may also be 
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applied to other forms of inter ferometric sensing. 
Therefore, the multi-wavelength illumination would occur 
simultaneously. To obtain the phase data there are a number 
of options the most common of which is phase stepping. Here 
5 the phase of each fringe pattern at each wavelength needs 
to be varied to give at least three image sets and hence 
some temporal acquisition is still required. The next most 
coiranon process is Fourier transform fringe analysis which 
derives the phase data from a single image. It will be 
10 appreciated that using the latter technique would enable 
all the required phase data to be captured in a genuinely 
''single snapshot" technique. 
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CLAIMS 

1. A multi-frequency interf erometry method for measuring the 
absolute phase of an electromagnetic wave, wherein the 

5 method comprises: 

selecting a target measurement range, L, within which 
absolute phase measurements are desired to be made, said 
phase measurements being related to a desired measurand of 
an object; 

10 determining a level of phase noise which will be 

present in wrapped phase measurements which will be made; 

for the selected target measurement range, and the 
determined phase noise level, calculating an optimum 
number, N, of wavelengths of electromagnetic radiation to 

15 be used in the multi-frequency interf erometry, where N^3; 

selecting an optimum series of values of said N 
wavelengths to achieve optimum noise immunity in wrapped 
phase measurements to be made; and 

carrying out multi-frequency interf erometry at the 

20 selected values of the N wavelengths so as to make at least 
one wrapped phase measurement at each of the N wavelengths, 
which wrapped phase measurements are processed to obtain an 
absolute phase measurement related to the desired measurand 
of the object. 

25 

2. The method according to claim 1, wherein the optimum 
number, of wavelengths is the minimum number of 
wavelengths required to obtain unambiguous phase 
measurements in the selected target measurement range, for 

30 the determined phase noise level. 

3. The method according to claim 1 or claim 2, wherein the 
method includes proposing a measure of process reliability 
associated with the determined phase noise level. 
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corresponding to a known probability of success in fringe 
order identification . 



4. The method according to claim 3, wherein the proposed 
5 process reliability is 6a, where a is the standard 
deviation noise in a discrete level heterodyne function 
generated from the wrapped phase measurements made at two 
different wavelengths . 

10 5. The method according to claim 3 or claim 4, wherein the 
proposed process reliability is used to calculate the 
minimum number, of wavelengths required to obtain 

unambiguous phase measurements in. the selected target 
measurement range . 

15 

6. The method according to any preceding claim, wherein the 
processing of the wrapped phase measurements includes 
heterodyne processing. 

20 7. The method according to claim 6, wherein the heterodyne 
processing produces a plurality of discrete level 
heterodyne functions each containing equal numbers of 
levels in the interval -ti to -^n. 



25 8 . The method according to any preceding claim, wherein the 
processing of the wrapped phase measurements includes 
iterative unwrapping . 



9. The method according to any of claims 1 to 5, wherein 
30 the selected values of said N wavelengths define a 
geometric series. 



10. The method according to any of claims 1 to 5, wherein 
the selected values of said N wavelengths are such that 
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they can be combined in a predetermined manner to generate 
a geometric series of synthetic wavelengths, 

11. The method according to claim 10, wherein the selected 
5 values of said N wavelengths are heterodyne processed to 

generate the geometric series of synthetic wavelengths. 

12. The method according to claim 11, wherein the selected 
values of the N wavelengths are chosen such that only one 

10 heterodyne operation is required to generate each synthetic 
wavelength in the geometric series. 

13. The method according to any of claims 10 to 12, wherein 
the values of said N wavelengths are selected in accordance 

15 with the following equation: 



where i=l,...,N, where N is the number of wavelengths, X.o is 

20 the wavelength associated with the largest frequency, and Xi 
is the wavelength of the ith frequency, and L is the target 
measurement range . 

14. The method according to any preceding claim, wherein 
25 the method is used for measuring the three-dimensional 
shape of an object which may have surface discontinuities, 
and the target measurement range is the range over which 
absolute depth measurements relating to the profile of the 
object are desired* 



15. The method according to claim 14, wherein the step of 



1/Xi = 1/A-o - (lAo) 



i-l/N-l 



(1/L) 



N-i/N-l 



30 
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inter ferometry 
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recording a series of fringe patterns obtained when 
the object is illuminated with electromagnetic radiation at 
each of the selected N optimum wavelengths; 

processing the recorded fringe patterns so as to 
5 obtain wrapped phase measurements in the form of a wrapped 
phase map for each of the N optimum wavelengths, and 
processing the wrapped phase maps to determine absolute 
fringe orders in the fringe patterns. 

10 16. The method according to claim 15, wherein the absolute 
fringe orders are calculated using at least one of 
heterodyne processing, recursive unwrapping and iterative 
unwrapping of the wrapped phase maps. 

15 17. The method according to claim 15 or claim 16, wherein 
at each said optimum wavelength a series of phase-stepped 
fringe patterns is recorded. 

18. The method according to any of claims 1 to 5, wherein 
20 the multi-frequency interf erometry is carried out using 
fringe projection, and the selected series of N optimum 
wavelengths are synthetic wavelengths projected on to the 
object. 

25 19. The method according to claim 18, wherein the selected 
values of the N synthetic wavelengths projected on to the 

object define a geometric series given by for 
/ = I,2,...,A , where X is number of wavelengths and N^^ is the 
number of projected fringes in the largest fringe set. 

30 

20. The method according to claim 18, wherein the values 
of the N synthetic wavelengths projected on to the object 
are selected in accordance with the following equation: 
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^f, =^/o -K)^. for /=i,2,...,^ 

where X is the number (X = N) of projected wavelengths, N^^ 
is the number of projected fringes in the largest fringe 
set/ and N^^ is the number of projected fringes in the i'* 
5 fringe set. 

21. The method according to claim 18, wherein the number, 
N, of wavelengths is equal to three, and the selected 
values of the three wavelengths are N^^, N^^-^N^ ^ ^nd ^N^^ -1, 

10 where N^^ is the number of projected fringes in the largest 
fringe set. 

22. The method according to any preceding claim, wherein 
the object is illuminated separately and sequentially with 

15 each of the selected optimum wavelengths of electromagnetic 
radiation. 

23. The method according to claim 14, wherein the object 
is illuminated with white light and a plurality of said 

20 fringe patterns are captured simultaneously by recording 
them with an image detector which can simultaneously 
capture image data at the selected optimum wavelengths. 

24. The method according to claim 23 wherein the image 
25 detector is a colour camera. 

25. The method according to any of claims 1 to 13, wherein 
the method is carried out in a single point ranging system 
in which the absolute phase measurements made are used to 

30 calculate absolute range. 
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26. The method according to claim 25, wherein the multi- 
frequency interf erometry is carried out by illuminating the 
object using a broadband femtosecond laser. 

5 27. The method according to claim 26, wherein the spectral 
range of the broadband pulses emitted by the laser is 
chosen so as to include the said optimum series of values 
of said N wavelengths. 

10 28. A method of measuring absolute fringe order in a phase 

measuring sensor, wherein the method comprises: 

selecting a number, N, of wavelengths of 

electromagnetic radiation to be used to illuminate an 

object/ where N>3; 
15 determining a level of phase noise which will be 

present in wrapped phase measurements which will be made; 
for the selected number, N, of wavelengths, and the 

determined phase noise level, selecting optimum values of 

the N wavelengths so as to achieve a maximum measurement 
20 range in which wrapped phase measurements can be made 

relating to a desired measurand of the object to be 

illuminated; 

recording a series of fringe patterns obtained when 
the object is illuminated with electromagnetic radiation at 
25 each of the selected N optimum wavelengths; and 

processing the recorded fringe patterns to obtain at 
least one wrapped phase measurement at each of the selected 
N optimum wavelengths, and processing the wrapped phase 
measurements to determine absolute fringe orders in the 
30 fringe patterns. 



29. A multi-frequency interferometer apparatus for shape 
measurement, comprising : 
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fringe projection means for generating a projected 
wavelength of illumination by projecting a predetermined 
fringe pattern on to an object; and 

image capture and recording means for capturing and 
5 recording a deformed fringe pattern obtained when said 
predetermined fringe pattern is projected on to the object^ 
the image capture means being disposed at an angle to the 
direction of illumination of the object with the projected 
fringe pattern; and 
10 data processing means for processing recorded deformed 

fringe patterns so as to obtain phase measurements 
therefrom; wherein 

the fringe projection apparatus is variable such that 
an optimum series of projected wavelengths may be 
15 generated/ the values of the wavelengths being such that, 
for a known level of phase noise in the apparatus, and a 
chosen number of projected wavelengths in said series, 
absolute fringe orders in the deformed fringe patterns are 
measurable over a maximum measurement range - 

20 

30. The apparatus according to claim 29, wherein the fringe 
projection means comprises a coherent fibre fringe 
projector for producing a pattern of Young's fringes across 
the object* 
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